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Abstract 

Tachyon-Brane inflationary universe model in the context of a Chaplygin gas equation of state 
is studied. General conditions for this model to be realizable are discussed. In the high-energy 
limit and by using an exponential potential we describe in great details the characteristic of this 
model. Recent observational data from the Wilkinson Microwave Anisotropy Probe experiment 
are employed to restrict the parameters of the model. 
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I. INTRODUCTION 



Inflationary universe models have solved many problems of the Standard Hot Big Bang 
scenario, for example, the flatness, the horizon, and the monopole problems, among others 
In addition, its has provided a causal interpretation of the origin of the observed 
anisotropy of the cosmic microwave background (CMB) radiation, and also the distribution 
of large scale structures. |3[]. 

In concern to higher dimensional theories, implications of string/M-theory to Friedmann- 
Robert son- Walker (FRW) cosmological models have recently attracted great deal of at- 
tention, in particular, those related to brane-antibrane configurations such as space-like 
branes[4|. The realization that we may live on a so-called brane embedded in a higher- 
dimensional Universe has significant implications to cosmology [5j. In this scenario the 
standard model particles are confined on the brane, while gravitations propagate in the 
bulk spacetimes. Since, the effect of the extra dimension induces additional terms in the 
Friedmann equation is modified at very high energies jf], Q], acquiring a quadratic term in 
the energy density. Such a term generally makes it easier to obtain inflation in the early 
Universe 0, For a review, see, e.g., Ref. 10]. 

In recent times a great amount of work has been invested in studying the inflationary 
model with a tachyon field. The tachyon field associated with unstable D-branes might be 
responsible for cosmological inflation in the early evolution of the universe, due to tachyon 
condensation near the top of the effective scalar potential 111 ], which could also add some 
new form of cosmological dark matter at late times Cosmological implications of this 
rolling tachyon were first studied by Gibbons 13] and in this context it is quite natural to 
consider scenarios in which inflation is drive by the rolling tachyon. 

On the other hand, the generalized Chaplygin gas has been proposed as an alternative 
model for describing the accelerating of the universe. The generalized Chaplygin gas is 



described by an exotic equation of state of the form p c h = —Ap c! ^, where p c h and 



Pch are 



the energy density and pressure of the generalized Chaplygin gas, respectively [14|. (3 is 



a constant that lies in the range < (5 < 1, and A is a positive constant. The original 
Chaplygin gas corresponds to the case j3 = 1 [la ]. Inserting this equation of state into the 



2 



relativistic energy conservation equation leads to an energy density given by [ijj] 



Pch 



A + 



B 



1+/9 



(1) 



G 3(l+/3) 

where a is the scale factor and B is a positive integration constant. 

The Chaplygin gas emerges as a effective fluid of a generalized d-brane in a (d+1, 1) 
space time, where the action can be written as a generalized Born-Infeld action \\M. These 

nnn 

models have been extensively studied in the literature [16|, H7|, 1181 ]. 

In the model of Chaplygin inspired inflation usually the scalar field, which drives inflation, 
is the standard inflaton field, where the energy density given by Eq.O^can be extrapolate for 
obtaining a successful inflation period with a Chaplygin gas model[19]. Recently, tachyon- 
Chaplygin inflationary universe model was considered in [20[ , and the dynamics of the early 
universe and the initial conditions for inflation in a model with radiation and a Chaplygin 
gas was studied in Ref . [2l| (see also |22J ) . 



A natural extension of Ref. 



201 ] is to consider the tachyon field as a degree of freedom 



on visible three dimensional brane. This work has been extended to include higher order 
corrections in slow-roll parameters and the formula has been widely used to confront this 
model with the observations. Moreover, we find constraints on the parameter A and the 
five-dimensional Planck mass or equivalently the brane tension. 

The outline of the paper is a follows. The next section presents a short review of the 
modified Friedmann equation by using a Chaplygin gas, and we present the tachyon-brane- 
Chaplygin inflationary model. Section II III deals with the calculations of cosmological per- 
turbations in general term. In Section HVl we use an exponential potential in the high-energy 
limit for obtaining explicit expression for the model. Finally, Sect JVl summarizes our findings. 



II. THE MODIFIED FRIEDMANN EQUATION AND THE TACHYON-BRANE- 
CHAPLYGIN INFLATIONARY PHASE. 

We consider the five-dimensional brane scenario, in which the Friedmann equation is 
modified from its usual form, in the following way 

where H = a/a denotes the Hubble parameter, represents the matter confined to the 
brane, k = 8ttG/3 = 8tt /3rrip (m p is the four-dimensional Planck mass), A 4 is the four- 
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dimensional cosmological constant and the final term represents the influence of bulk gravi- 
tons on the brane, where £ is an integration constant (this term appears as a form of dark 
radiation). The brane tension A relates the four and five- dimensional Planck masses via 
m p = a/3M|/(47tA), and is constrained by the requirement of successful nucleosynthesis as 
A > (lMeV) 4 23|]. We assume that the four-dimensional cosmological constant is set to 
zero, and once inflation begins the final term will rapidly become unimportant, leaving us 



with 



101 ] . Hence, the modified Friedmann equation reads 



H 2 



K 



A + 



,(1+0) 



1 

1+/9 



A + 



,(1+j8) 



1+/3 



2A 



(3) 



Here, p^ becomes p^ = V ((/))/ y 1 — 2 , and V{4>) = V is the tachyonic potential. Note that, 
in the low energy regime [A + p^ 1+/3 - ) ] 1 /( 1 +^) ^ A, the tachyon-Chaplygin inflationary model 
is recovered [20j, and in a very hight-energy regime, the contribution from the matter in 
Eq.Q becomes proportional to [A + p( 1+ ^] 2 /( 1+ #l in the effective energy density. 

We assume that the tachyon field is confined to the brane, so that its field equation has 
the form 

rh ■ V 

0, (4) 



+3H<j) + — 
l-(j) 2 V 

where dots mean derivatives with respect to the cosmological time and V = dV(4>)/d(f>. For 
convenience we will use units in which c — H — 1. 

The modification of the Eq.Q is realized from an extrapolation of Eq.([T]), where 
the density matter p m ~ a -3 in introduced in such a way that we may write p c h = 



A + p 



m 



1+/3 



A + 



1+/3 



and thus, we identifying p m with the contributions 
of the scalar tachyon field which gives Eq.([3]). The generalized Chaplygin gas model ma y b e 
viewed as a modification of gravity, as described in Ref.j^, for chaotic inflation, in Ref. 19]. 
and for tachyon- Chapligin inflationary universe model in the low-energy limit, in Ref. 20]. 
Different modifications of gravity have been proposed in the last few years, and there has 
been a lot of interest in the construction of early universe scenarios in higher-dimensional 
models motivated by string/M-theory 25) . It is well-known that these modifications can 
lead to important changes in the early universe. In the following we will take (3 = 1 for 
simplicity, which means the usual Chaplygin gas. 

During the inflationary epoch the energy density associated to the tachyon field is of the 



4 



order of the potentia 
and <j> < 3H(p 



13 



i.e. 



V. Assuming the set of slow-roll conditions, i.e. 2 1 



261 ]. the Friedmann equation (j3J) reduces to 
H 2 ^kVA + V 2 



1 + 



+ v 2 



2A 



and Eq. (j3J) becomes 



V 



(5) 



(6) 



Introducing the dimensionless slow-roll parameters 23], we write 



# m 2 



H 2 167T 



/2 



1 + 



(A+V 2 ) 1 / 2 



V 2A 



and 
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m l 












m 2 


2HV ~ 


16 7T 



1/" 



{ ^ (A + V 2 ) 1 / 2 



V 



I 2 



1 + 



2A 

(A + V 2 ) 1 / 2 
2A 



(7) 



(8) 



(9) 



20] 



Note that in the low-energy limit, yA + p 2 <C A, the slow-parameters are recovered 
The condition under which inflation takes place can be summarized with the parameter 
e satisfying the inequality e < 1, which is analogue to the requirement that a > 0. This 
condition could be written in terms of the tachyon potential and its derivative V, which 
becomes 



V' 2 



1 | (A + V 2 ) 1 / 2 
A 



< ^1 {A + V 2 f 2 
mi 



1 + 



(A + V 2 ) 1 ' 2 



-, 2 



2A 



(10) 



Inflation ends when the universe heats up at a time when e ~ 1, which implies 



(A + V 2 ) 1 / 2 

1 + ^ 



167T 



mi 



(A + Vff 2 



1 + 



(A + V 2 ) 1 / 2 



2A 



(11) 



However, in the high-energy limit [A + 



-211/2 



[A + V 2 } 1 / 2 > A Eq.JIIJ becomes 



y/2 ^4vr (A + V 2 ) 2 



A 



The number of e-folds at the end of inflation is given by 



P ^0 



vM+172 
2A 



(12) 
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or equivalently 



N = — 



8tt f y f V VA + V 2 



m 



y/2 



1 + 



VA + V 2 
2A 



dV. (13) 
-(4n/m 2 p \)tff[V (A + 



Note that in the high-energy limit Eq. (1131) becomes N ~ 
V 2 )/V l2 ]dV. 

In the following, the subscripts * and / are used to denote the epoch when the cosmo 
logical scales exit the horizon and the end of inflation, respectively. 



III. PERTURBATIONS 



In this section we will study the scalar and tensor perturbations for our model. It was 
shown in Ref. 28j that the conservation of the curvature perturbation, 1Z, holds for adiabatic 
perturbations irrespective of the form of gravitational equations by considering the local 
conservation of the energy-momentum tensor. However, we note here that even though 
the effect of bulk to the cosmological perturbations can not be trivially negligible, it can 
be shown that the main correction of the spectrum in the brane-world inflation is just the 



modification of the slow-roll parameters 
spectrum of the curvature perturbations is given Vn — ( ^-r 



291 ] (see also [30(]). For a tachyon field the power 



where Z s 



V(l-< 



-3/2 



271 ] . Under the slow- roll approximation, the power spectrum of curvature perturbations is 



estimated to be 



27j 



H 2 V 1 



128tt (V(A + V 2 y' 2 



1 + 



(A + V 2 ) 1 / 2 
2A 



2n(hJ V 3////. V V' 2 
Note that in the low-energy limit the amplitude of scalar perturbation given by Eq.([T 
coincides with Ref. 20]. 



(14) 



The scalar spectral index n s is given by n s — 1 = ^n^rr, where the interval in wave number 
is related to the number of e-folds by the relation d\nk((f)) = —dN(<f>). From Eq.( fl4l) . we 
get, n s ~ 1 — 2(2e + 7 — 77), or equivalently 



m 



4ir 
V' 2 



V -(A + V 2 )- 1 ' 2 



1 + 



A 



(A + V 2 ) 



, (A+y 2 ) 1 / 2 

1 ~r 2A 



(A + v 2 y/ 2 

2A 



V' 2 V" 
2V 2 ~~ V 



-, -1 



(15) 
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One of the interesting features of the five-year data set from Wilkinson Microwave 
Anisotropy Probe (WMAP) is that it hints at a significant running in the scalar spec- 
tral index dn s /dhik — a s [3|. From Eq. (ll5j) we get that the running of the scalar spectral 
index becomes 



4{A + V 2 ) 

W' 



1 + 



(A+V 2 ) 1 / 2 



2A 



1 + 



(A+V 2 ) 1 / 2 



[2e,* + 7,*-f7,*] £ - 



(16) 



In models with only scalar fluctuations the marginalized value for the derivative of the 
spectral index is approximately —0.03 from WMAP-five year data only [9(. 

On the other hand, the generation of tensor perturbations during inflation would produce 
gravitational waves and this perturbations in cosmology are more involved since gravitons 



propagate in the bulk. The amplitude of tensor perturbations was evaluated in Refs. 
and 1 321 



3l| 
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V g = 2Ak (JL\ 2 F 2 (x) ^^n 2 (A + V 2 ) 1 / 2 



7T Z 



1 



(A + V 2 ) 1 ' 2 



2A 



F 2 (x) 



(17) 



where x = Hm p ^3/ (4-7rA) and 



F{x) — Vl + x 2 — x 2 sinh l (l/x) 



-1/2 



Here the function F(x) appeared from the normalization of a zero-mode. The spectral index 
n is given by n a = -r-r-S- = — at ~it — w- 

9 o J 9 d\nk ff^i ^l+x 2 

From expressions (TT^I) and fflTj) we write the tensor-scalar ratio as 



r{k) 



' 9 

Pn 



k=k„ 



V' 2 F 2 (V) 



(18) 



k=k t 



.3 k V (A + V 2 ) [1 + (A + V 2 y/ 2 /2X} 2 
Here, is referred to k = Ha, the value when the universe scale crosses the Hubble horizon 
during inflation. 

Combining WMAP five-year dataQ with the Sloan Digital Sky Survey (SDSS) large 



scale structure surveys [33j, it is found an upper bound for r given by r(k* ~ 0.002 Mpc 



-v 



< 



0.28 (95%CL), where k* ~0.002 Mpc 1 corresponds to / = r k ~ 30, with the distance to the 
decoupling surface tq= 14,400 Mpc. The SDSS measures galaxy distributions at red-shifts 
a ~ 0.1 and probes k in the range 0.016 h Mpc _1 < k <0.011 h Mpc" 1 . The recent WMAP 
five-year results give the values for the scalar curvature spectrum P-ji(k^) ~ 2.4 x 10~ 9 and 
the scalar-tensor ratio r(k*) = 0.055. We will make use of these values to set constrains on 
the parameters appearing in our model. 



IV. EXPONENTIAL POTENTIAL IN THE HIGH-ENERGY LIMIT. 



Let us consider a tachyonic effective potential V(4>), with the properties satisfying 
V{4>) — ► as — > oo. The exact form of the potential is V(<p) = (1 + acp) exp(— a0), 
which in the case when a — > 0, we may use the asymptotic exponential expression. This 
form for the potential is derived from string theory calculations ll|, ]34j. Therefore, we simple 
use 

V{<f>) = V e- a< f>, (19) 

where a and Vq are free parameters. In the following we will restrict ourselves to the case 
in which a > 0. Note that a represents the tachyon mass 26|, |35(. In Ref. 12) is given 



an estimation of these parameters in the low-energy limit and A — ■> 0. Here, it was found 
Vq ~ I0~ 10 rrip and a ~ lCT 6 m p . In the following, we develop models in the high-energy 



limit, i.e. y/A + V 2 > A. 

From Eq. fTl3|) the number of e- folds results in 

4tt 



N 



X a 2 m 2 



[KV*)-h{V f )l 



where 



V 2 

h(V) = — + A lnV. 



(20) 



(21) 



On the other hand, we may establish that the end of inflation is governed by the condition 
e — 1, from which we get that the square of the tachyonic potential becomes 



v{4 = ttf 



V = — 

s 8vr 



A a 2 m\ - 8 n A + X a 2 m 2 p (X a 2 m 2 p - 1Q ti A) , (22) 



and 



<f) f 



a m Tl 



Note that in the limit A 



2Vf 

we obtain Vf 




3tt 
am, 



v / A/(2 v /tt) and 
coincides with that reported in Ref. 12] . 

From Eq. (114p we obtain that the scalar power spectrum is given by 



(23) 

l/y/3, which 



Pn(k) 



16vr 



3 m® a 2 A 3 



(A + V 
V 



2\3 



and from Eq. (|18|) the tensor-scalar ratio becomes 



r(k) 



4 m 2 X 2 a 2 

7T 



V 

{A + V 



2\2 



F 2 {V) 



(24) 



(25) 
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By using, that V = —a V, we obtain from Eq.(jT5 



m; 

1 ~ - 



A a 2 



4tt (A + V 2 ) l(A + V 2 ) 



AV 2 



and from Eq. ffTBT) that 



\ 2 a 4 m 4 r 



4vr 2 



3V 2 -5A 
(A + V 2 ) 4 



V 2 . 



(26) 



(27) 



The Eqs. (1241) and ( 1261) has roots that can be solved analytically for the parameters a and 
A, as a function of n s , P-^., V and A. The real root solution for m 2 , and A becomes 



a 



2vr 



and 



where 



4 4 I/ 5 + 6 2 P n (n s -l)V 2 \ 2 m 4 + X(3P n (l- n s ) A 2 m 4 - 8 2 V 3 ) 



A=^(2^ 2 -K) 



(28) 



(29) 



K = ^16^ 4 + 3(n s - l)P n V\ 2 m 4 . 

From Eq. (f29|) and since A > 0, the ratio V 3 /A 2 satisfies the inequality y 3 /A 2 < (1 — 
ns)P-jim 4 /4:. This inequality allows us to obtain an upper limit for the ratio V^ 3 (0)/A 2 
evaluate when the cosmological scales exit the horizon, i.e. K* 3 /A 2 < 2.4 x 10~ u m 4 . Here, 
we have used the WMAP five year data where P-ji(k^) ~ 2.4 x 1CT 9 and n s (k*) ~ 0.96. 

One again, note that in the limit A — > 0, the constrains a ~ 7 x 10~ 3 M 5 and V* ~ 



4 x 10 _4 M| are recovered 12]. Here, we used the relation m p = M| a/3/(47tA). 

In Fig. 1 we have plotted the adimensional quantity A 2 /m^ versus the adimensional scalar 
tachyon potential evaluated when the cosmological scales exist the horizon V^/m 4 . In doing 
this, we using Eq.( |271) that has roos that can be solved for the brane tension A, as a function 
of a s , m, A and V. For a real root solution for A, and from Eqs. ( 1281) and ( 1291) we obtain a 
relation of the form A = /(K) for a fixed values of a s , n s and Pu- In this plot we using the 
WMAP five year data where P n (K) ^ 2.4 x KT 9 , n s (K) ~ 0.96, a s (K) ~ -0.03. In Fig. 2 
we have plotted the tensor-scalar ratio given by Eq. fT25l) versus the adimentional parameter 
A/m,p. The WMAP five-year data favors the tensor-scalar ratio r ~ 0.055 and the from 
Fig. 2 we obtain that A parameter becomes A ~ 2.6 x 10~ 25 mp. In Fig. 3 we have plotted 
the tensor-scalar ratio given by Eq. (|25l) versus the adimentional parameter a 2 /m 2 . We note 
that for r ~ 0.055 we obtain a 2 ~ 1.3 x 10~ 12 m 2 . 
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For these values of the A and a parameters we get the values V* ~ 1.3 x 10 



-12™, 4 



ni 



V f ~ 8.9 x 10~ 14 m^ and A ~ 5.1 x lO^m 4 ~ 4 x 10~ 5 M 5 4 . Also, the number of e-folds, N, 
becomes of the order of N ~ 52.7. We should note also that the A parameter becomes smaller 
by two order of magnitude and the a parameter becomes similar when it are compared with 
the case of tachyon-Chaplygin inflation in the low-energy limit [3]. 



x10' 




V./m 



x10 



FIG. 1: The plot shows the adimentional square of the brane tension (A/m 4 ) versus the 
adimentional scalar potential V*/rrip. Here, we have used the WMAP five-year data where 
Pn(h) ^ 2.4 x 1(T 9 , n s (h) ~ 0.96 and a s {h) ~ -0.03. 



V. CONCLUSIONS 



In this work, we have studied the tachyon-Chaplygin inflationary model in the context of a 
branewold scenario. In the slow-roll approximation we have found a general relation between 
the scalar potential and its derivative. This has led us to a general criterium for inflation to 
occur (see Eq.t flOi) ). We have also obtained explicit expressions for the corresponding scalar 
spectrum index n s and its running a s . 
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FIG. 2: The plot shows the tensor-scalar ratio r versus the adimentional parameter A/rrip. Here, we 
have used the WMAP five-year data where Pn(K) 2.4xl0~ 9 , n s {K) ~ 0.96 and a s (h) ~ -0.03. 



By using an exponential potential in the high-energy regime and from the WMAP five 
year data, we found the constraints of the parameters A and a from the tensor-scalar 
ratio r (see Figs. 2 and 3). In order to bring some explicit results we have taken the 
constraints A ~ 10~ 25 mp and a ~ 10 _6 m p , from which we get the values K ~ lCr 12 mp, 
Vf ~ lCF 13 mp, A ~ 10 _13 mp and ~ 53. Here, we have used the WMAP five year data 
where P n {K) ^ 2.4 x 1CT 9 , n s (K) ~ 0.96, a s {K) ~ -0.03 and r(K) ~ 0.055. Note that the 
restrictions imposed by currents observational data allowed us to establish a small range for 
the parameters that appear in the tachyon-brane-Chaplygin inflationary model. 

We have not addressed reheating and transition to standard cosmology in our model (see 
e.g., Ref. |36i]). Specifically, it will be very interesting to know how the reheating temperature 
in the hight-energy scenario, contributes to establish some constrains on the parameters of 
the model. We hope to return to this point in the near future. 
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FIG. 3: The plot shows the tensor-scalar ratio r versus the adimentional parameter a 2 /nip. Here, 
we have used the WMAP five-year data where Pn{k*) — 2.4 x 1CT 9 , n s (k*) ~ 0.96 and a s (k*) ~ 
-0.03. 

Acknowledgments 

This work was supported by the "Programa Bicentenario de Ciencia y Tecnologia" 
through the Grant N° PSD/06. 



[1] A. Guth, Phys. Rev. D 23, 347 (1981). 

[2] A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett. 48, 1220 (1982); A complete description of 
inflationary scenarios can be found in the book by A. Linde , Particle physics and inflationary 
cosmology (Gordon and Breach, New York, 1990). 

[3] J. Dunkley et al. [WMAP Collaboration], arXiv:0803.0586l [astro-ph]; G. Hinshaw et al, 
larXiv:0803.0732l [astro-ph]. 

[4] A. Sen, JHEP 0204, 048 (2002). 



12 



[5] K. Akama, Lect. Notes Phys. 176, 267 (1982); V. A. Rubakov and M. E. Shaposhnikov, Phys. 
Lett. B 159, 22 (1985); N. Arkani- Hamed, S. Dimopoulos, and G. Dvali, Phys. Lett. B 429, 
263 (1998); M. Gogberashvili, Europhys. Lett. 49, 396 (2000); L. Randall and R. Sundrum, 
Phys. Rev. Lett. 83, 3370 (1999); 83, 4690 (1999). 

[6] P. Binetruy, C. Deffayet, and D. Langlois, Nucl. Phys. B565, 269 (2000) ; P. Binetruy, C. 
Deffayet, U. Ellwanger, and D. Langlois, Phys. Lett. B 477, 285 (2000). 

[7] T. Shiromizu, K. Maeda, and M. Sasaki, Phys. Rev. D 62, 024012 (2000). 

[8] R. Maartens, D. Wands, B. A. Bassett, and I. P. C. Heard, Phys. Rev. D 62, 041301 (2000). 

[9] J. M. Cline, C. Grojean, and G. Servant, Phys. Rev. Lett. 83, 4245 (1999); C. Csaki, M. 
Graesser, C. Kolda, and J. Terning, Phys. Lett. B 462, 34 (1999); D. Ida, JHEP 0009, 014 
(2000); R. N. Mohapatra, A. Perez-Lorenzana, and C. A. de S. Pires, Phys. Rev. D 62, 105030 
(2000); R. N. Mohapatra, A. Perez-Lorenzana, and C. A. de S. Pires, Int. J. Mod. Phys. A 



16, 1431 (2001); Y. Gong, arXiv:gr-qc/0005075 



[10] R. Maartens, arXiv:gr-qc/0101059 



[11] A. Sen, Mod. Phys. Lett. A 17, 1797 (2002). 

[12] M. Sami , P. Chingangbam and T. Qureshi, Phys. Rev. D 66, 043530 (2002). 
[13] G. W. Gibbons, Phys. Lett. B 537, 1 (2002). 

[14] M. C. Bento, O. Bertolami and A. Sen, Phys. Rev. D 66, 043507 (2002). 
[15] A. Kamenshchik, U. Moschella and V. Pasquier, Phys. Lett. B 511, 265 (2001); N. Bilic, 
G. B. Tupper and R. D. Viollier, Phys. Lett. B 535, 17 (2002). 



[16] H. B. Benaoum, [arXiv:hep- th/0205140; A. Dev, J.S. Alcaniz and D. Jain, Phys. Rev. D 67 
023515 (2003); G. M. Kremer, Gen. Rel. Grav. 35, 1459 (2003); R. Bean, O. Dore, Phys. 
Rev. D 68, 023515 (2003); Z.H. Zhu, Astron. Astrophys. 423, 421 (2004); H. Sandvik, M. 
Tegmark, M. Zaldarriaga and I. Waga, Phys. Rev. D 69, 123524 (2004) ; L. Amendola, I. 
Waga and F. Finelli, JCAP 0511, 009 (2005); P.F. Gonzalez-Diaz, Phys. Lett. B 562, 1 
(2003) ; L.P. Chimento, Phys. Rev. D 69, 123517 (2004); L.P. Chimento and R. Lazkoz, Phys. 
Lett. B 615, 146 (2005); U. Debnath, A. Banerjee and S. Chakraborty, Class. Quant. Grav. 
21, 5609 (2004); W. Zimdahl and J.C. Fabris, Class. Quant. Grav. 22, 4311 (2005); P. Wu 
and H. Yu,Class. Quant. Grav. 24, 4661 (2007); M. R. Setare, Phys. Lett. B 654, 1 (2007); 
M. Bouhmadi-Lopez and R. Lazkoz, Phys. Lett. B 654, 51 (2007). 

[17] M.C. Bento, O. Bertolami and A. Sen, Phys. Lett. B 575, 172 (2003); M.C. Bento, O. Berto- 



13 



lami and A. A. Sen, Phys. Rev. D67, 063003 (2003). 
[18] M. Makler, S.Q. de Oliveira and I. Waga, Phys. Lett. B 555, 1 (2003); J.C. Fabris, S.V.B. 

Goncalves and P.E. de Souza, |arXiv:astro-p h/0207430; M. Biesiada, W. Godlowski and M. 

Szydlowski, Astrophys. J. 622, 28 (2005); Y. Gong and C.K. Duan, Mon. Not. Roy. Astron. 

Soc. 352, 847 (2004); Y. Gong, JCAP 0503, 007 (2005) . 
[19] O. Bertolami and V. Duvvuri, Phys. Lett. B 640, 121 (2006). 
[20] S. del Campo and R. Herrera, Phys. Lett. B 660, 282 (2008). 
[21] G. A. Monerat et al, Phys. Rev. D 76, 024017 (2007). 

[22] R. Herrera, Phys. Lett. B 664, 149 (2008); S. del Campo and R. Herrera, Phys. Lett. B 665, 
100 (2008). 

[23] J. M. Cline, C. Grojean and G. Servant, Phys. Rev. Lett. 83, 4245 (1999) 
[24] T. Barreiro, A. A. Sen, Phys. Rev. D 70, 124013 (2004). 

[25] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999); T. Shiromizu, K. Maeda and M. 
Sasaki, Phys. Rev. D 62, 024012 (2000); D. Dvali, G. Gabadadze and M. Porrati, Phys.Lett. 
B 485, 208 (2000); K. Freese and M. Lewis, Phys. Lett. B 540, 1 (2002); R. Maartens, Lect. 
Notes Phys. 653 213 (2004); A. Lue, Phys. Rept. 423, 1 (2006). 

[26] M. Fairbairn and M. H. G. Tytgat, Phys. Lett. B 546, 1 (2002). 

[27] J. C. Hwang and H. Noh, Phys. Rev. D 66, 084009 (2002). 

[28] S. Tsujikawa, D. Parkinson and B. A. Bassett, Phys. Rev. D 67, 083516 (2003). 

[29] K. Koyama, D. Langlois, R. Maartens and D. Wands, JCAP 0411, 002 (2004). 

[30] K. Koyama, A. Mennim and D. Wands, Phys. Rev. D 72, 064001 (2005); K. Koyama, 
S. Mizuno and D. Wands JCAP 0508, 009 (2005); A. Cardoso, K. Koyama, A. Mennim, 
S. S. Seahra and D. Wands, Phys. Rev. D 75, 084002 (2007); K. Koyama, A. Mennim, 
V. A. Rubakov, D. Wands and T. Hiramatsu, JCAP 0704 001 (2007); K. Koyama, A. Men- 
nim and D. Wands, Phys. Rev. D 77, 021501 (2008). 

[31] D. Langlois, R. Maartens and D. Wands, Phys. Lett. B 489, 259 (2000). 

[32] M. Sami, N. Savchenko and A. Toporensky, Phys. Rev. D 70, 123528 (2004) 

[33] M. Tegmark et al., Phys. Rev. D 69, 103501 (2004). 

[34] D. Kutasov, M. Marino and G. W. Moore, JHEP 0010, 045 (2000). 

[35] A. de la Macorra, U. Filobello and G. German, Phys. Lett. B 635, 355 (2006) 

[36] B. Feng and M. Li, Phys. Lett. B 564, 169 (2003); C. Campuzano, S. del Campo, and R. Her- 

14 



rera, Phys. Rev. D 72, 083515 (2005), [Erratum-ibid. D 72, 109902 (2005)]; C. Campuzano, 
S. del Campo and R. Herrera, Phys. Lett. B 633, 149 (2006); C. Campuzano, S. del Campo 
and R. Herrera, JCAP 0606, 017 (2006); S. del Campo and R. Herrera, Phys. Rev. D 76, 
103503 (2007); J. C. Bueno Sanchez and K. Dimopoulos, JCAP 11, 007 (2007). 



15 



